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Abstract. 

We study the effects of optically thick clumps, non-void inter-clump medium, 
variation of the onset of clumping, and velocity dispersion inside clumps on the forma- 
tion of resonance lines. For this purpose we developed a full 3-D Monte Carlo Radiative 
Transfer (MCRT) code that is able to handle 3-D shapes of clumps and arbitrary 3-D 
velocity fields. The method we developed allows us to take into account contributions 
from density and velocity wind inhomogeneities to the total opacity very precisely. The 
first comparison with observation shows that 3-D density and velocity wind inhomo- 
geneities have a very strong influence on the resonance line formation, and that they 
have to be accounted for in order to obtain reliable mass-loss rate determinations. 



1. Introduction 

There is a lot of observational and theoretical evidence that the winds of hot stars are 
inhomogenoues. Stocha stic variable stru ctures in the He u A 4686A emission line in 



^ Pup were revealed bv lEversberg et all (.19 98). and explained as an excess emission 
from the wind clumps. iMarkova et al.l ( 20051) investigated the Ho- line-profile variability 



for a large sample of 0-type supergiants. They concluded that the properties of this 
variability can be explained by a wind model consisting of clumps. 

Understanding of clumping i s a prerequisite to obtain correct mass-loss rates from 
modeling the stellar spectra (see iFuUerton et all l2006l : iBouret et al.l l2005l : IPuIs et al" 



12006 : Oskinova et al. 2007) . 

A possible way how to study wind clumping is to use simplifications. The ap- 
proximate treatment of clumping is usually based on the "microclumping" (clumps are 
optically thin at all frequencies) or "macroclumping" (clumps may be optically thick 
for certain wavelengths) approaches. Common additional assumptions are the void in- 
terclumped medium (ICM) and the monotonic velocity field. The enhancement of the 
density inside clumps is usually describe d by the s o called "clumping factor" D. 

In the most recent studies of Sundq vist et al.l ((2010, 2011) it was shown that the 
detailed density structure, non-monotonic velocity field, and the ICM are all important 
for the line formation in the massive hot star winds. In order to be able to study the 3-D 
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nature of wind clumping and to take into account all properties of wind clumping men- 
tioned here, we developed a realistic full 3-D MCRT code (Surlan et al., in preparation). 
Here we present the basic effects of the different clump properties on the resonance line 
formation and present the first preliminary comparison with observation. 



2. The clumped wind model 

Our wind model is grid-less and does not require any symmetry. Instead of using a 
predefined grid, we introduce an adaptive integration step Ar. First we create a snapshot 
of clumps inside the wind and then we follow photons along their paths using the Monte 
Carlo approach. The photon frequencies are expressed in the local co-moving frame. 
All distances are expressed in units of the stellar radius (/?*). 

The velocity field can be arbitrary in our wind model. However, here we assumed 
a radial monotonic velocity v,- field and the standard yS-velocity law for the underlying 
smooth wind {Vr = vp). The velocity inside clumps ii^. is expressed as 

r - 

Vcir) = v/s(r.) - Vdisir) — — ^, (1) 

where v^hir) = m vp{r) describes the velocity dispersion using the deviation parameter 
m (0 < m < 1). The monotonic velocity at the position is vpir':), and Z,- and r': are the 

radius and the absolute position of the center of the /-th clump, respectively. 

The opacity of the wind is calculated using the parameterization given by lHamami 

(Il980h . 



-Yur <Px, <Px = —p. 



X{r) ^ -Yv£ <Px\ (px^—e- , (2) 



where xo is the opacity parameter that corresponds to the line strength, Vco is the termi- 
nal velocity of the wind, and (p^ is the absorption profile, which depends on the dimen- 
sionless frequency x = (v/vq - l)c/vo, where vo is the Doppler-broadening velocity 
which comprises both the thermal broadening and microturbulence. The ionisation is 
assumed to be constant. 

For simplicity we assume a spherical shape of the clumps with a radius that varies 
with the distance r from the star, / = /(r). A clump radius is determined using the clump 
separation parameter Lq as 



l(r) = Iq yfrhvip); where 1q = Lq^ (3) 

w{r) = v,./vco is the velocity in the units of the terminal speed iioo, and D is the clumping 
factor. The average clump separation L is also depth dependent and given as L(r) = 
Lo-y'r^ w{r). The density inside clumps is assumed to be by the clumping factor D 
higher than the smooth wind density. In the first approximation, D is assumed to be 
independent on r and given by 

L3(r) 

D = , , D > 1. (4) 
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The ICM density is determined by the interclump density parameter d which is 
also assumed to be depth independent and < d < 1. If the factor d = 0, the ICM 
is void, otherwise the ICM is rarefied by the factor d compared to the smooth wind 
density. For the case of dense clumps and void ICM, all the mass of the wind is put 
into clumps, and the volume filling factor is fy = l/D. For the case of dense clumps 
and non-void ICM, the mass of the wind is distributed between clumps and ICM and it 
follows that /y - (1 - d)/{D - 1). 

The total number of the clumps is determined by the condition of the mass con- 
servation inside the wind. Clumps are created one by one and when the total volume 
of created clumps (Vd) is larger than the total volume of the wind (Vw) multiplied by 
fv, i.e. if Vci > fvVw, then the number of created clumps represents the total number 
of clumps in our snapshot clump distribution. 



3. Effect of different clump properties 

By varying the model parameters we examined how macroclumping (defined by Lq 
and D), non mono tonic velocity field (vdis), non void ICM (d), and onset of clumping 
(rci) may influence weak (j^o = 0.5), intermediate (xo = 50), and strong (^q = 500) 
resonance line profiles. Calculations are performed assuming ~ 10"* clumps inside 
wind, which corresponds to Lq = 0.5 for the wind with r = 10. The clumping factor 
is set to D = 10. The Doppler-broadening velocity is vd = 50 km s~^ A snapshot of 
distribution of the clumps is shown in Fig.[T] upper right panel. 

The first set of line profiles (CI in Fig. [T|l is calculated assuming that clumping 
starts at the surface of the star (r^ = 1) with void ICM (d = 0), and monotonic velocity 
field (iidis - t^/j)- The weak (upper left panel), intermediate (lower left panel), and 
strong (lower right panel) lines show strong line strength reduction compared to the 
smooth wind (full lines). If the wind is clumped, more photons can escape through 
"holes" between clumps. This leads to lower effective opacities, and, consequently, 
less absorption. The line strength reduction depends on the number of clumps in the 
wind. The absorption dip close to Vco is caused by lower probability of photon escape 
through lower number of velocity "holes". 

In order to investigate how the onset of clumping influences line profiles, we cal- 
culated the second set of line profiles (C2 in Fig. [B assuming that clumping is only 
above = 1.3. Other model parameters are the same as for the first set of line profiles 
calculations. A very strong and broad absorption near the line centers appears. It is due 
to a higher absorption of the smooth part of the wind, which is in the region 1 < r < r^. 
Such deep absorption is not observed and it disappears if clumping starts closer to the 
wind base (perhaps even from the surface of the star) or if the ICM is sufficiently dense. 
The weak line almost reproduces the smooth wind model in this case. For a small value 
of xo, the clumps are optically thin and they have a very small influence on the line 
profile compared to the smooth wind. 

Varying the ICM density parameter d, a significant impact on the line profiles was 
found, especially on strong lines. Keeping the same parameters as for the second set of 
line profile calculation, but choosing d = 0.05 we calculated the third set of line profiles 
(C3 in Fig.[T]l. While the weak line is unaffected, the intermediate and strong ones show 
an increase of absorption in the blue edge of the line and a decrease of absorption near 
the centre of the line. The degree of the line saturation depends on the ICM density. For 
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Figure 1 . Line profiles for weak (upper left panel), medium (lower left panel), and 
strong (lower right panel) lines for the clump distribution in the upper right panel. S 
- the smooth model; Cn denotes the clumped models - CI (rd = 1), C2 (rd = 1.3), 
C3 (rci = 1.3, c/ = 0.05), C4 (r^ = 1.3, = 0.05, Udis = 0.2 v/i). 

a particular value of d, saturation of strong lines and desaturation of the intermediate 
lines can be achieved. 

Taking into account the non-monotonic velocity field, an absorption at velocities 
higher than Vco appears. Setting y^is - 0.2 v/^ and keeping other model parameters the 
same as for the third set of line profiles, we calculated the fourth set of line profiles (C4 
in Fig. [T]l. The higher i;dis> the larger is the velocity span of clumps, and, consequently, 
there are more velocity overlaps (the "holes" in the velocity field are smaller) and the 
probability of the photon escape is lower. 



4. Comparison with observations 

As the first attempt to compare with observations, we try to fit the P v doublet of ^ Pup 
observed by the copernicus satellite (thin line in Fig. O. As the terminal velocity of 
^ Pup we adopt Vco = 2250 km s~\ and the value offi = 0.9. The synthetic spectrum of 
the smooth and clumped wind models are calculated. The model parameters are chosen 
to fit the observed spectrum best. It can be seen that the predicted unclumped (smooth) 
P-Cygni profile (full line with crosses) of P v is much stronger than the observed one 
(thin full line). However, the synthetic spectrum for the clumped model (thick full line) 
fits the strength of the observed line very well. Therefore, using the unclumped model 
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Figure 2. Comparison of line profiles calculated using different wind models with 
COPERNICUS observations (thin full line). The unclumped wind model is shown by the 
full line with crosses, while clumped model is shown by the thick full line. 



can lead to underestimating the empirical mass-loss ra tes. T hese results are consistent 
with .Oskinova et al..(2007.) and .Sundqvist et al..(,20ial201lh . 



5. Conclusions 

Using our 3-D Monte Carlo Radiative Transfer code, we investigate the influence of 
the 3-D stellar wind clumping on the resonance line formation. Our results show that 
clumping may significantly reduce the strength of lines and, if not included, may lead 
to a wrong mass-loss rate determination. The 3-D wind nature, the onset of clumping, 
non void ICM, and non-monotonic velocity field are all important model ingredients. 
Our models allow to obtain these parameters from the comparison between model and 
observed lines, and thus gain deep insight in the structure of stellar winds, and obtain 
reliable mass-loss rate measurements. 
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